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Abstract 

We study the role of a direct intracavity polarization coupling in the dynam- 
ics of transverse pattern formation in type-II optical parametric oscillators. 
Transverse intensity patterns are predicted from a stability analysis, numeri- 
cally observed, and described in terms of amplitude equations. Standing wave 
intensity patterns for the two polarization components of the field arise from 
the nonlinear competition between two concentric rings of unstable modes in 
the far field. Close to threshold a wavelength is selected leading to stand- 
ing waves with the same wavelength for the two polarization components. 
Far from threshold the competition stabilizes patterns in which two different 
wavelengths coexist. 
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I. INTRODUCTION 



Pattern formation is an ubiquitous manifestation of nonlinearity Jiy^] which presents 
specially interesting features in nonlinear optical systems |3|-|5]]. The search of transverse 
structures in nonlinear optical systems is actively pursued for several reasons which include 
their possible application in all-optical signal processing and the investigation of macroscopic 
manifestations of quantum phenomena. These structures are the result of the interaction of 
nonlinearity and diffraction in transverse spatially extended devices such as nonlinear optical 
cavities of large Fresnel number. Among the nonlinear systems analyzed, optical parametric 
oscillators (OPO's) have received a lot of attention from the theoretical viewpoint. Avail 
able results include the analysis of pattern formation [jSj-TTOlj, noise sustained structures [11 



domain walls |T2}-|T4|] and localized structures [|To]-|T^|] . A growing interest in these transverse 
structures in OPOs arise also from the study of quantum spatial correlations present in these 
patterns |5|,|TP|-|2"T| . Transverse patterns in OPO have been recently observed |22[| . 

In an OPO two first harmonic (FH) fields (signal and idler) are generated inside the 
crystal by parametric down conversion of the external pump field. In type-I OPO signal and 
idler fields have the same state of linear polarization. In type-II OPO they are orthogonally 
polarized. This polarization degree of freedom can be used for a nonlinear construction of 
new states of the emitted light. For example, by means of a direct polarization coupling 
produced by an intracavity quarter wave-plate (A/4-plate) it is possible to produce p3lj2~4}| 
states in which signal and idler are degenerate in frequency and phase locked. A general 
question that we address in this paper is the effect of this type of direct polarization coupling 
in the problem of transverse pattern formation in type-II OPO. 

When considering transverse spatial degrees of freedom in type-II OPO without a direct 
polarization coupling there are two different regimes. In one of them, characterized by a 
positive effective detuning homogeneous solutions are selected at threshold with an arbi- 
trary relative phase between signal and idler fields. For effective negative detuning a finite 
wave number is selected at threshold and a phase pattern -traveling wave (TW)- is asymp- 



totically selected for each FH field, while the intensity remains homogeneous p5|,26|. The 



effects of direct polarization coupling between signal and idler in type-II OPO for a positive 
effective detuning was discussed in refs. [^,^8|: Spatial domains of equivalent, but different, 
self-phase locked homogeneous solutions appear. They are separated by phase polarization 
domain walls. These walls are generally Bloch walls whose motion leads to complex spa- 
tiotemporal states. In this paper we consider the situation of negative effective detuning in 
order to characterize the way in which polarization coupling modifies the process of pattern 
formation. In particular one expects that the coupling between signal and idler can generate 
standing waves (SWs) -i.e, stripe intensity patterns- from the TWs that exist for signal and 
idler when there is no direct polarization coupling. 

For a type-II OPO with an intracavity A/4-plate and for FH's negative detunings we 
predict, and numerically confirm, that there is a threshold of pattern formation above which 
the FH's far fields exhibit during the transient dynamics two concentric rings of growing 
unstable modes. Therefore, the coupling does not only lead from TWs to SWs, but it 
introduces two different wavelengths, giving rise to an interesting problem of wavelength 
competition and pattern selection. The values of the different wavelengths are controlled by 
the A/4-plate polarization coupling. In a "symmetric case" (i.e, all the dynamical parameters 
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for signal and idler are equal) we show that the two wavelengths coexist for long times. 
Real and imaginary parts of each FH field present in this case SW patterns with different 
wavelength which consist in domains of stripe patterns. For asymmetric FH's coefficients 
and near threshold, we show that one of the wavelengths dominates and intensity stripe 
patterns -standing waves- with the same wavelength emerge for signal and idler. Far from 
threshold the dynamics is similar to the symmetric case. This transition from wavelength 
coexistence to dominance of one of the two competing wavelengths is described by means 
of an amplitude equation analysis. The amplitude equations give a full description of the 
observed dynamics. 

The paper is organized as follows. Section 2 reviews the mean field equations for this 
system. In section 3 we analyze the threshold for pattern formation. The instability is 
characterized in terms of the eigenf unctions associated with the critical modes. In section 4 
we discuss numerically the dynamics of transverse pattern. In section 5 we derive amplitude 
equations for the critical modes, which give a clear interpretation of the observed dynamics. 
We summarize our main conclusions in section 6. 



II. MEAN FIELD EQUATIONS 

We consider an optical parametric oscillator that consists of a ring optical cavity filled 
with a birefringent, nonlinear quadratic medium and externally pumped by a uniform laser 
beam. A direct polarization coupling between the FH fields, that take into account the effect 
of an intracavity quarter wave-plate is also included in the model. This wave plate provides 
a polarization mixing between the signal A x and idler A y fields p3|^4| . The signal and 
the idler can be either frequency degenerate or non-degenerate, depending on the frequency 
selection rules imposed by the combined effects of the parametric down-conversion, the cavity 
resonances and phase- matching [^-31|, but they are always polarization non-degenerate 



(type-II interaction). In the mean field approximation, and considering the paraxial and the 
single longitudinal mode approximation for all the fields, the equations describing the time 
evolution for the linear polarization components of the second harmonic (B X)y (x,y,t)) (SH) 
and the first harmonic (A Xty (x,y,t)) (FH) slowly varying envelopes of the electric fields, in 
a type-II, phase-matched OPO are [[23],[^J : 

d t B x = i x [-{\ + iA' x )B x + ia' x V 2 B x + 2iK A x A y + E + c'B y ] 

d t B y = tJ[-(1 + iA' y )B y + ia' y V 2 B y - c'*B x ] 

d t A x = 7x[-(l + iA x )A x + ia x V 2 A x + iK A*B x + cA y ] 

d t A y = 7l/ [-(l + iA y )A y + ia y V 2 A y + iK A* x B x - c*A x ] (1) 

where with no loss of generality we take A x , B x as ordinary polarized beams and A y , B y 
as extraordinary polarized f2"6|j . The coefficients J x ,y,l xy ( cav hy decay rates), A X)yi A' xy 



(cavity detunings) and a x>yi a' x (diffraction coefficients) are defined as in refs. P,p7fl; due 



to the birefringence of the nonlinear crystal all the diffraction coefficients can be slightly 
different, even when the signal and idler are frequency degenerate. Other parameters are 
the nonlinearity K and the injected pump E (bifurcation parameter) that, for the sake of 
simplicity, we take it real and polarized along the same direction than the phase-matched 
component of the second harmonic field B x . Hence, B y neither is pumped nor is nonlinearly 
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coupled with other components of the field. Under these conditions B y does not influence the 
dynamics although it is linearly coupled to B x . In eqs.(|l|) the direct polarization coupling 
constants (c, d) account for the effects produced by the A/4-plate and they are related to 
the phase mismatch and the axes of the A/4-plate by: 

c ~ sin(20) exp(if) (2) 

where (f) is the angle between the A/4-plate's fast axis and the principal axis of the crystal. 
The phase £ is the round-trip phase shift between the signal and idler at frequency degeneracy 
(or between B x and B y for d). The coupling strength |c| depends on a number of factors, 
including mode matching and Poynting vector walk-off . Here we assume perfect phase- 
matching and propagation along the optical axis (i.e. we neglect spatial walk-off). We focus 
on the effects produced by the rotation angle which is an important experimental parameter 
to control the effects described below. We note that other forms of similar linear coupling 
terms between signal and idler considered in |2"7| , |2"8| are associated with a birefringent and/or 
dichroic cavity mirror in type-II OPO. 

The linear coupling c breaks the phase invariance, that eqs.(|l|) have for c = 0, under 
changes of the relative phase of the FH fields: [A x , A y ] — ► [exp(icp) A x , exp(-iip) A y ]. How- 
ever, the phase of the signal and idler fields can be adjusted in order to include the phase 
of c in the FH fields; i.e, given c = |c| exp(i£), the transformation: [Az, A/] — > [A^Ad = 
[exp(— i£/2) A x , exp(i£/2) Ay] leaves eqs.([l|) unchanged except for the replacement c — > \c\. 
For the sake of clarity we will then present our main results for c real. The generalization 
for c complex is trivial using the phase transformation of the FH fields just described. 



III. LINEAR STABILITY ANALYSIS 

A. Threshold analysis 

In this section we present the linear stability analysis of the steady-state solution of 
eqs.(P corresponding to the OPO operating below the threshold of signal generation. This 
trivial uniform steady state (off-state) is: 

A X = Ay = 

b x = (i + iA' y ) e /[i - a^a; - |cf + i(A4 + a;)] 

B y = -d*B x /(l + iA' y ) (3) 

The threshold for transverse pattern formation is determined by linearizing eqs.(^) around 
this solution and looking for instabilities. The steady state becomes unstable only along the 
directions of the FH components (A x , A y ) of the eigenvectors, and thus the analysis reduces 
to the study of two linearly-coupled complex equations. Because of the complex nature of 
the field variables, it is convenient to consider the real and the imaginary parts of these 
equations for each FH field. The most general solution of the perturbations is hence given 
by a linear superposition of terms of the form: [Re(A Xty ), Im(A Xjy )] ~ exp[iqr + X(q)t] where 
X(q) is the growth rate of the perturbations and q is its transverse wave vector. 

For c = d = 0, the linear stability analysis shows that the trivial solution is stable for 
\F\ < \F C \, where F is a normalized pump intensity 



4 



F = K E Q /{l + iA' x ) 



(4) 



For A = 7^ A x + j y Ay > the most unstable mode corresponds to an homogeneous solution 
go = 0. In this paper we focus on the case A < for which the unstable modes at threshold 
(|F C | = 1) correspond to transverse traveling waves A X ,A* ~ exp[iqf + X(q)t], whose two- 
dimensional wave-vector q lies on a circle centered at with radius go = \l — A/a, where 
a = •jxCCx + ^yOLy For |F| > |F C | and A < any g mode on the circle, and the opposite 
mode for the orthogonal component of the field, can be selected at threshold by means of 
spontaneous symmetry breaking. Hence a phase pattern appears above threshold for A x and 
A y -traveling waves-, with opposite wavevector, while the intensity remains homogeneous in 



both polarizations [26,25]. It should be noted that when idler and signal fields are degenerate 
both in frequency and polarization -i.e, type-I DOPO-, eqs.(|IJ) must be solved with the 
further condition A x = A y . In this case standing waves states are selected at threshold ||. 
The linear coupling of A x with A y considered in this paper is expected to produce standing 
waves for each polarization component also in type-II OPO. 

For the threshold for pattern formation remains at \F C \ = 1, but a main difference is 
that now the wavevectors of the most unstable modes lie on two concentric circles of different 
radius. This gives rise to a wavenumber competition in the process of pattern formation, as 
shown below. Closed expressions are hard to obtain in the general case, but for 'jx = ly = 1 
the two eigenvalues of the linearized equations that characterize the instability are 



Ai, 2 = 7 



where we have defined: 



-1 + - J A \F\ 2 - 4 |c| 2 - 2 (62 + 92) ± 2T 



(5) 



Qj = Aj + acj q 2 (j = x, y) 
T : 



'-4 |F|2 (Q x - Q y y + 4 [el 2 (9, + 9,)2 + (92 - 92)2 
and we have introduced the normalized pump amplitude: 

F = (1 + iA' y ) K E /[l - A' x A' y - |cf + i(A' x + A' y )) 



(6) 



(7) 



which coincides with F given by eq.(|4]) for d = 0. Two other eigenvalues remain always 
negative. 

To avoid cumbersome expressions, analytical results are derived in this section for the 



particular case ^ x = 7 y = 7, a x 
eigenvalues Ai^g) become: 



a and A^, = A^ = A (< 0). In this case, the 
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•1 + 



(A + a|g-1 2 ± c ) 



where the plus (minus) sign corresponds to Ai (A2). From eq. 
instability for perturbations with an arbitrary wavevector q: 



A + a |gf ± c 



(8) 

we get the threshold of 



(9) 



where the plus (minus) sign corresponds to F\ (F2). Therefore, for c < —A the instability 
takes place at the critical threshold |F C | = 1 and the unstable modes at threshold correspond 
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to transverse traveling waves whose two-dimensional real wave- vector q lies on either of two 
concentric circles, centered at with radius q\ 2' 



9? >2 = 9i,2 = — (10) 

In figure 1 we show the instability threshold for perturbations of different wavenumbers. 
The threshold F c = 1 is the same for c = and c 7^ 0, but for c 7^ the instability takes 
place at two different wavevectors of modulus qx j2 indicated in the figure. Homogeneous 
perturbations (q = 0) have a larger instability threshold. The homogeneous phase locked 
solutions associated with this threshold are discussed in the appendix. We will not consider 
here the case c > — A > for which at F c — 1 only the mode q 2 becomes unstable. In this 
case the mode q± = becomes unstable for larger values of the pump. 

The values of the most unstable wavenumbers q\ and q2 depend on the absolute value of 
c as follows from eq. fllPp (figure 2). In particular, q\ 2 coincide with q\ for c = and vary 
linearly with |c|. The polarization coupling splits the circle of unstable modes for c = in 
two circles. The value of c controls the magnitude of the split. We remark that for c ~ |A|, 
|(fi| ~ and patterns with very large wavelength can be expected. For \F\ > \F C \ there is a 
band of unstable modes associated with each eigenvalue. In figure 3 we show the real part of 
the eigenvalues Ai,2 as a function of |g| for the critical case (\F\ = \F C \ = 1) and for one case 
above threshold. All the modes with i?e(Ai ) 2) > are linearly unstable. The wavevectors 
with modulus \q\\ or \q 2 \ have the same maximum growth rate: Ai 2(51,2) = 7 ( — 1 + 1-^1)- 

When the damping, detuning or diffraction parameters for signal and idler are different 
the competing modes of wavenumbers q\ and q 2 have different thresholds. It follows from 
the numerical analysis of eq.(f|) that the smaller wavenumber q\ becomes first unstable at 
\F C \ — 1. In figure 4a we show the growth rate for perturbations of the trivial state as a 
function of \q\ for a case in which \F\ = 1.0019. The unstable mode of smaller wavenumber 
qi dominates in the linear regime. The difference in growth rates for q\ and q 2 depend on 
the pumping level. For larger values of the pump, the growth rates of both unstable modes 
are of the same order, being the growth rate of q\ larger than the one of q 2 (see fig. 4.b). 
In the next section we show that this fact deeply affects the nonlinear mode competition 
dynamics of the system. 

We finally note that for asymmetric FH's coefficients, the critical modes have in general 
real eigenvalues (i.e. Im(Xx,2 ((7*1,2)) = 0). However, for very small values of c, i.e., when 
|<fi| ~ \<h\ there is a small interval of values of the external pump F, which includes the 
critical value F c = 1, where Im(Xi )2 (qi^ 2 )) 7^ 0. In this limit of c — > the instability becomes 



convective, similarly to a situation considered in ref. p2 . 



B. Critical modes 

Next we consider some features of the early time dynamics of pattern formation that can 
be understood in terms of the eigenvectors corresponding to the eigenvalues of the linearized 
problem discussed above. First, we introduce the far field as the Fourier transform of the 
near field, where the near field is the transverse field configuration at the input/output cavity 
mirror. The far field components A^{t) of A x (for example) are defined by: 
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1 f°° ~ 

A x(r,t) = — Ajj(t) exp(zgr) dq x dq y (11) 

Z7T J-oo 



In figure 5 we show numerical results |33| for a typical transverse profile of the A x and 
A y FH fields at an early time after the pump is increased beyond its threshold value. These 
results correspond to the case of symmetric coefficients in which two competing wavenumbers 
have the same growth rate. In Fig. 5a we show the near field of the signal intensity pattern 
and its far field. The two concentric rings of the far field correspond to unstable wavectors q 
with arbitrary orientation and wavenumber around q\ (inner ring) and q2 (outer ring). The 
near field is the result of the interference among all the unstable modes of both rings in the 
far field. However, the interference takes place in such a way that the real part of A x is 
associated with the unstable modes of the outer ring, while its imaginary part is associated 
with the unstable modes of the inner ring. In addition there is a high correlation between 
the transverse structures observed in the signal A x and idler A y fields. This is illustrated in 
Figs. 5b, 5c and 5d. In Figs. 5b and 5c we show the near field of the real and imaginary 
parts of A x and A y . It is observed that the real and imaginary parts of the field A x support 
transverse patterns with different wavelength. The same fact is observed in the idler field 
A y , but real and imaginary part have a different wavelength than for A x . In fact we observe 
that Re(A x ) ~ Im(A y ), while Im(A x ) ~ Re(A y ). In Fig. 5d we show the far field of Re(A x ) 
and Im(A x ). This gives evidence of the different wavenumber associated with Re(A x ) and 
Im(A x ). The two rings in the far fields of Re(A x ) and Im(A x ) correspond, respectively, to 
the outer and inner rings of Fig. 5a. 

These numerical facts can be explained in terms of the eigenvectors A 12 (^) associated 
with the eigenvalues Ai j2 introduced in eq.(H). They can be written as 

[Re(A x ), Im(A x ), Re(A y ), Im(A y )f = At(q) = C\ [ki, 1, 1, -«if exp(iqf) 

A 2 (g) = C 2 [1, -«2, «2, if exp(iqr) (12) 



where we have introduced: 



«i 



— Ai + F - I 
A + a \q\ 2 + c 



-X2 + F-I 

k 2 = — — 13 

A + a \q\ 2 — c 



and the normalization constants Cj = J 2 (1 + |/%| 2 ), j = 1,2. For simplicity we take 
here F to be real (i.e. resonant pump field). Exactly at threshold, K\ 7 2 vanish for the 
corresponding critical wave- vector: Ki >2 (91,2) = 0. The general dependence of on the 
wavenumber at threshold is shown in Fig. 6. At threshold F = F c = 1 and the eigenvectors 
^■1,2(0) are damped for any q except Ax(<f = q\) and A 2 (g = q 2 ): 

A 1 (q 1 ) = -L [0,1,1, Of exp(i^f) 

A 2 (g 2 ) = i= [1, 0, 0, if exp(ig 2 r) (14) 

which are marginal (zero growth rate) and define the direction in the functional space along 
which the instability takes place. The form of the eigenvectors (0) explain our numerical 
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finding (Fig. 5) that at short times, when the linear approximation to the dynamics remains 
valid, the components of the FH's fields (Re(A x ), Im(A y )) and (Im(A x ), Re(A y )) only sus- 
tain patterns with wavevector of modulus |^| or \q\\ respectively. As we discuss in the next 
section, this gives rise to pattern formation with competing wavelengths that can be rather 
different. 

There is an interesting symmetry in the far field components of the first harmonics 
fields. We focus in the A x field, but the discussion is also valid for A y . Our numerical 
results indicate, as shown for example in Fig. 7, that the real (imaginary) part of the far 
field component A^(t) is an odd function for q — q\ and an even function for q = q 2 (even 
function for q = q± and an odd function for q = q 2 ). This property implies that 

\ (15) 

Given that in the linear regime at threshold only excitations with wavectors q 2 or q\ con- 
tribute to the resulting structures in the real- valued fields Re(A x ) or Im(A x ), respectively, 
we find that the basic excitations for Re(A x ) and Im(A x ) are standing waves of the form 

Re(A x ) ~ Aq 2 exp(iq 2 r) + A_ g - 2 exp(-iq 2 f) 
Um(A x ) ~ Aft exp(wfif) + A-^ exp(-igif) (16) 

These standing waves are the interference between two opposite modes of the same ring 
of the far field which satisfy eq. (|15D . Both modes have the same amplitude but the global 
phase of the superposition of the two modes is different in each ring. 

It is interesting to note that the circularly polarized components A± of the vectorial FH 
field give a natural description of the instability. We have already discussed that, at threshold 
Re(A x ) ~ Im(Ay), and Im(A x ) ~ Re(A y ), so that A + = (A x + iA y )/^/2 = iy/2 Im(A x ) and 
A_ = (Ax — iA y )/y/2 = y/2 Re(A x ). Therefore, it follows from eq.fllBD that the instability 
for A + (A_) takes place at qi (q 2 ). The two circularly polarized components will emerge 
at the instability as standing waves intensity patterns of different wavenumber. In fact, for 
symmetric coefficients and in the linearized version of eqs. (1) around the trivial solution 
(0), A + and A_ are decoupled. They are nonlinearly coupled through the pump field B x . 

Our above discussion is for real values of the parameter c. When c is complex (£ ^ 0) 
both the real and imaginary part of A x and A y have contributions of unstable modes of 
wavenumber q± and q 2 . Therefore we observe transverse patterns with competing wave- 
lengths in the real and the imaginary parts of A XjV . However, by changing the global phase 
of the FH fields, the problem can be considered in terms of c real, as previously pointed out. 
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IV. TRANSVERSE INTENSITY PATTERNS 



In this section we give a numerical description |33] of the patterns that are asymptotically 
selected after a regime of nonlinear competition among the unstable modes of wavenumber 
qi and q 2 . A theoretical justification of these results is given in the next section in terms of 
an amplitude equation. We consider separately the cases of symmetric and nonsymmetric 
coefficients for the FH fields. 



S 



A. Symmetric coefficients 



When j x = 7^ = 7, a x = a y = a and A x = A y = A the growth rates of the most unstable 
modes in the two circles of radius qi and q 2 are equal and the instability takes place at both 
circles simultaneously (see Fig. 3). In this case the nonlinear competition keeps wavevectors 
in both circles in the far field excited for long times. The real and imaginary part of the 
FH fields show patterns with different wavelengths, as discussed in section 3. To illustrate 
the dynamical evolution we show in Figs. 8 and 9 snapshots of the pattern configuration at 
two different times. The two rings in the far field persist at long times. The pattern which 
appears in the long time dynamics presents domains in which real and imaginary parts of 
each of the FH near fields show standing waves of arbitrary orientation and of different 
wavelength for the real and imaginary parts. The arbitrary orientation of these standing 
waves comes from the spontaneous choice of two opposite wavevectors in the corresponding 
ring of the far field. Thus, the general picture is that there is long time competition among 
standing waves of different orientations and two different wavenumbers. 

A limiting situation is the case in which the inner ring collapse to the point q — 0. This 
situation takes place when c = \A X \. In this case, the growth of uniform domains occurs in 
one of the components of the FH vector field while the other component sustains standing 
waves of local arbitrary orientation as shown in figure 10. This structure represents a case in 
which patterns and uniform domains coexist in the same complex field as an effect induced 
by the direct polarization coupling. This situation takes place near signal resonance for very 
small values of c. 

B. Nonsymmetric coefficients 

When the damping, diffraction or detuning coefficients of signal and idler are different, 
the nonlinear mode competition depends very much on how far above threshold the OPO is 
pumped. Near threshold the relative difference in the magnitude of the growth rate of the 
unstable modes on the circles of wavenumbers qi and qi is important, as follows from Fig. 
4. a. This fact produces a strong change in the dynamics of the system. In Figs. 11, 12 and 
13 we show snapshots of configurations at different times of the dynamical evolution. After 
a transient linear regime, discussed in section 3 and which is represented here by Fig. 11, 
there is a nonlinear competition between the two circles of unstable modes at intermediate 
times. This is shown in Fig. 12 where the far field of Re(A x ) and Im(A y ) are seen to 
have competing contributions from the two circles. At late times the inner circle wins the 
competition and the final pattern is a standing wave of wavenumber qi both for A x and A y . 
This nonlinear wavenumber selection can be traced back to the behavior of the growth rate, 
as shown in the next section. For long times there is also a spontaneous breaking of the 
rotational symmetry and a standing wave in an arbitrary direction is selected, as shown in 
Fig. 13. The real and imaginary part of the FH fields show in this case patterns with the same 
wavelength. The resulting structures originate in the interference between two arbitrary -but 
opposite- wavevectors of the inner circle. The resulting stripe intensity pattern is similar 
to that predicted for type-I DOPO in the sense that it is the interference between two 
opposite traveling waves ||. However, physically, the energy and momentum conservation 
in the parametric down conversion of pump photons only implies off-axis emission of idler 
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and signal photons along two opposite directions without interference between them because 
they have orthogonal polarization. Due to polarization coupling, pure traveling wave are 
not solutions of eqs. ([!]). The A/4-plate provides a mixing of polarization that allows that 
photons of the same FH field interfere producing a standing wave. This phenomenon is 



the same that occurs in resonantly coupled complex Ginzburg-Landau equations f55fl : the 
linear (polarization) coupling allows the formation of standing waves as the result of the 
interaction between opposite modes in the far field -in this case, the inner ring-. 

Far from threshold, the growth rates of the unstable modes qi and q 2 are of the same 
order (see Fig. 4.b) and the dynamics of the system in this regime is equivalent to the 
symmetric case, as we prove in the next section. Therefore, in the asymmetric case the 
intensity of the external pump can be used to stabilize both rings of the far field which 
have a competing coexistence for long times and far from threshold. In figure 14 we show a 
typical long time state for this regime. 



V. AMPLITUDE EQUATIONS 

Close to the instability threshold, and using the general methods of nonlinear dynamics 
and pattern formation theory one may derive amplitude equations for the patterns described 
numerically in the previous section. As it is now well known, these equations are able 
to describe pattern evolution, selection and stability. In particular, if, for simplicity, one 
considers only the critical modes identified in section 3, field variables may by expressed as: 

[Re(A x ), Im(A x ), Re{A y ), Im(A y )f = A 1 [0, 1, 1, 0] T exp(zg 1 f) + A_i [0, 1, 1, 0] T exp(-igi 
+ A 2 [1, 0, 0, l] r exp(iq 2 f) + A_ 2 [1, 0, 0, 1] T exp(-iq 2 r) (1 

where Ai(t) are slowly evolving amplitudes associated with the modes with wavevector 
<fj. Given the relations flT5| ) we have chosen here to define the amplitudes A±i, A± 2 as 
A±< ri = iA ±1 and A±$ 2 = A± 2 

Standard analysis 0] leads to the following evolution equations for the symmetric case 
(Q x = Q y and d = A' x = for simplicity) 

d t A ±l = 7 (|F| - l)At! + ^-[A ±1 (B x + B* x ) + A Tl (B x + B*)± 2?1 + %A 2 (B x - B*) ±q ^ ?2 

+ iA_ 2 (B x - B*)±j 1+ f 2 ] 
d t A ±2 = 7 (|F| - l)A ±2 + ^[A ±2 (B x + B* x ) + A t2 (B x + B*) ±2?2 + iA x (B x - B* x ) ± ^ 

+ l A_ 1 {B x -B x ) ± ^ 2 ] (18) 

where 

d t (B x ) ? = -(1 + ia'q 2 ) (B x ) € + 2iK {A^A y ) ? (19) 

a' = a' x and (~.)q indicates Fourier transform. Using eqs.(|l^) and time scales separation 
near threshold (\F\ -1 < 1) one may write {B x )j = Bj as an expansion in the amplitudes 
Ai. One has, at lowest order: 
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B = -AK (AiA-i + A 2 A„ 2 ) 

2K A± i A± i 
±2q ^~ l + Aia'qf 
£±a±& * (20) 

This leads to: 

d t A ±1 = 7 (|F| - 1)A ±1 - AK 2 A 2 ±l A Tl - AK 2 A ±1 A 2 A„ 2 - A^A^A 2 ^ 

d t A ±2 = 7 (|F| - 1)A ±2 - AK 2 A 2 ±2 A T2 - AK^A^A^ - Ar, 2 KlA T2 A 2 ±2 (21) 

where 77^2 = 2 [1 + (4a' g 2 2 ) 2 ] 

Since A_i = A\, A^ 2 = A 2 , as follows from (|15|) and the definition of the amplitudes 
A±\-, A± 2 after eq. (Ti7|), one may finally write: 

c^! = 7 (|F| - 1)A, - 4K 2 A 1 [(1 + r/OIAxI 2 + \A 2 \ 2 } 

d t A 2 = 7 (|F| - 1)A 2 - AKl A 2 [(1 + m )\A 2 \ 2 + |^| 2 ] (22) 

These equations describe a weak competition between the amplitudes A\ and A 2 , so that 
the stable steady states of the system are given by: 



l^i 



2 V2 7(1*1-1) 



?7i + % + ^7i^72 4if, 
r/i 7(1^ 



^| 2 = ' 2 " (23) 



f/i + ^2 + 771772 4K, 



These expressions are still valid when A' and d are different from zero, with 77^2 given by: 

= 1 x 1 + c' 2 + (A- + 4a' qj,) 2 

Vl ' 2 2 [l + (c' + A' + 4a'g 2 i2 ) 2 ][l + (c'-A'-4a'g 2 2 ) 2 ] 1 ' 



The asymptotic states of the dynamics, eqs. (|23[ ), correspond to patterns built on two 
wavevectors q\ and q 2 , as confirmed by the numerical results of the previous section. 
However, when the kinetic coefficients of x and y field components are slightly different 
(Q x 7^ Q y , \Q X — Q y \ <C Q x + ©y) the growth rates Ai j2 of critical modes become different 
(see fig. 4) and the amplitude equations (^) become: 

d t A, = X 1 A 1 - AKl A, (IA| 2 + |^ 2 | 2 ) - 477^2 A x \A,\ 2 

3 t A 2 = \ 2 A 2 -AKlA 2 {\A 1 \ 2 + \A 2 \ 2 )-Ai l2 KlA 2 \A 2 \ 2 (25) 

where Ai > A 2 . In this case, patterns with wavevector q\ only develop at, and slightly 
beyond threshold. The corresponding steady state is 

w = 4iWTri> w = (26) 

Such patterns are stable provided A 2 < Ai/(1 +771). When Ai reaches A 2 (l +771) they become 
unstable and in the domain where Ai is larger than A 2 (l+77i) patterns with both wavevectors 
qi and q 2 are stable. Their amplitude is given by: 
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Ai (l + tfc) - A 2 




4^(771 + 772 + r?ir? 2 ) 
A 2 (I + 771) - Ax 



(27) 



4^(771 + 7/2 + 7?i77 2 ) 



The two types of behavior (close and far from threshold) have been described in our 
numerical results in the previous section. Furthermore, on increasing pumping beyond 
threshold, one observes a crossover between monomode patterns, with wavevectors cor- 
responding to the maximum growth rate, and bimodal ones. For growth rates given 
by eqs.(|5|), transition between mono and bimodal patterns occurs, for small anisotropies 



In conclusion we have shown that standing wave intensity patterns can be generated in 
type-II optical parametric oscillators. They appear spontaneously in the transverse plane 
when there is a direct polarization coupling between the signal and idler fields, produced 
for example by an intracavity quarter wave-plate. Such a coupling gives also rise to two 
competing wavelengths in the system. 

In the transient dynamical regime after the pump is switched-on above its threshold value, 
there is a competition between two rings of unstable modes. This gives rise to transverse 
patterns with different wavelengths for the real and imaginary part of the FH's fields. We 
have described two dynamical regimes. In the first one, which correspond to symmetric 
FH's parameters, the far field is composed by two concentric rings. Real and imaginary part 
of each of the FH's fields shows patterns with different wavelength. For asymmetric FH's 
coefficients, the dynamical regime depends on the distance to threshold. Near threshold, the 
dynamical process of pattern competition leads to the dominance of a unique wavelength, 
selecting a transverse stripe intensity pattern with the same wavelength for signal and idler. 
Far from threshold, the dynamics is equivalent to the symmetric case, so that the strength 
of the external pump can be used to stabilize the two competing wavelengths. Amplitude 
equations for the vectorial (critical) modes have been derived and they confirm our numerical 
observations. In particular, it is worth noting that the structure of these equations reflect 
the vectorial nature of the fields, which introduces nontrivial couplings between the modes. 
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VIII. APPENDIX: UNIFORM PHASE-LOCKED SOLUTIONS 



In addition to the trivial uniform solution given by eq.(£J), eqs.(|IJ) have other uniform 
stationary solutions for A x /A y > considered in ref. p4j. These are the dominant solu- 



tions when transverse effects are not taken into account. For the sake of completeness, 
a brief description of these solutions for A XjV < is given here. Eqs. ([!]) admit for 
A x ,y = a Xt y exp(i(p Xj y) two uniform stationary solutions, which take the form: 



cos(4> d + £) 



cos 



i-r 2 
2| c |r 

(A x + TA y ) - 2 |c| T sin(0 d + O 



2 K E c p 

2 - <',, /'-() ' 



2c„£ n r sm(<j>i) - 1 - r 2 



4c p r 2 



a y = Ta x (29) 

where T 2 = A x /A y , c p = 1/(1 — |c'| 2 ) and 4>d,s = 4>x -F 4>y The existence of these solutions 
requires for r ^ 1 that |c| > |1 — r 2 |/(2T). This relation defines a circle in the complex 
plane of c inside which no stationary uniform solutions exist. 

The first of eqs.(^) indicates that these homogeneous solutions are self phase-locked. 
Due to the fact that the function arccos is a multi- valued function in the range [— tt, 7r] , two 
branches of uniform phase- locked solutions exist. Each branch has a different threshold. 
Uniform solutions are never spontaneously observed for A x ^ y < because they have a larger 
threshold of instability than the solution with a finite wavenumber. The threshold of the 
uniform branch with lower threshold is indicated by a diamond in figure 1. Note that this 
is larger than the threshold for pattern formation \F C \ = 1. The second branch of uniform 
solutions, with a larger threshold value, is unstable and it is not observed, not even for 
positive FH's detunings where uniform solutions dominate and domain walls between them 
have been reported P7| , P8] . The homogeneous solutions described here are different of those 
in Fig. 10 obtained in the limit of q\ = 0. 
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FIGURES 
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FIG. 1. Scaled threshold of instability F for the trivial stationary solution (eq.Q) as a function 
of q 2 for the case of symmetric FH's coefficients. Solid (dashed) curve gives the threshold value of 
F for c / (c = 0). The critical threshold \F C \ = 1 is indicated as a dotted line. The instability 
takes place at two different wavectors indicated by q\ and q2 (qo denotes the unstable wavector 



for c = 0). Parameters are Kq = 1, A a 



-0.8, A' = 



^0 — J-j ^x,y — — U.O, L± Xjy — U, a x>y — ^. x y 

and d = 0.01. The diamond indicates the threshold of instability for homogeneous perturbations 
(q = 0). 
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\q\ 

FIG. 3. Real part of the eigenvalues \\^ as a function of \q\ for the case of symmetric FH's 
coefficients. The solid lines correspond to F = 1.06 and the dotted to criticality, F = F c = 1. Left 
branches correspond to Ai and the right ones to A2. Here 7^ = j' xy = 1 and the values of the 
other parameters are the same as in Fig. 1. The most unstable modes q\ and qi correspond to the 
maximum of each line. The level -Re(Ai 5 2) = is also indicated as a dashed horizontal line as a 
reference. 
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FIG. 4. Growth rate of the unstable modes as a function of \q\ in a case in which the FH's 
coefficients are different. Fig. 4. a) Near threshold: F = 1.0019 (solid line) and F = 1 (dotted line). 
Fig. 4.b) Far from threshold: F = 1.02. Coefficients values are Kq = 1, A X)J , = —0.8, A' x = 0, 
lx = 7 ' x = 0.9901, ly = i y = 1.01, a x = 0.2475, a y = 0.2525, a' XjV = 0.125, c = 0.4 and d = 0.01. 
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Figure 4.b 
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FIG. 5. A snapshot at t = 900 of the FH's fields spontaneously generated from random initial 
conditions close to the trivial steady-sate given by (3). a) left: near field \ A X \; right: far field \A^\. 
b) A x field. Left: Re(A x ); right: Im(A x ). c) A y field. Left: Im(A y ); right: Re{A y ). d) Absolute 
value of the Fourier transform. Left: Re(A x ); right: Im(A x ). Parameters are F = 1.16, Kq = 1, 
A x>y = -0.8, A' x>y = 0, lx , y = i Xty = 1, a x , y = 0.25, a' x<y = 0.125, c = 0.4 and d = 0.01. 
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FIG. 6. Amplitudes (solid line) as a function of \q\. The real part of the critical eigenvalues 
Ai,2 (dotted lines) and the zero level (dashed line) are also plotted as reference. Here F = 1 and 
the values of the other parameters are the same of figure 3. 
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FIG. 7. A typical result for a cut of the far field (FF) amplitude A^{t) along the line q y = for 
the A x field during the transient state, a) Real part; b) Imaginary part. Parameters are F = 1.06, 
K = 1, A x>y = -0.8, A' x>y = 0, 7x = 1 ' x = 0.9901, ly = i y = 1.01, a x = 0.2475, a x = 0.2525, 
' - 0.125, c = 0.4 and d = 0.01. 
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Figure 7.b 



22 



FIG. 8. Snapshots as in Fig. 5 at time t = 7800. Parameters are the same of Fig.5 




23 




24 





FIG. 10. A snapshot at time t = 3000 of the field A x spontaneously generated from random 
initial conditions for a case in which the mode q = is one of the most unstable modes. Figures a), 
b), c) and d) show respectively the intensity, far field, real part and imaginary part of the signal 
field. In c) and d) patterns with very different wavelength can be appreciated. Parameter values 
are F = 1.08, K = 1, A x>y = -0.6, A! xy = 0,j x = i x = ~f y = j' y = 1, a x = a y = 0.25, c = 0.6 and 
d = 0.01. 
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FIG. 11. Snapshots as in Fig. 5 at time t = 100. Parameter values are Eq = 1.002, Kq = 1, 
A X:V = -0.8, A' xy = 0, lx = j' x = 0.9901, ~/ y = i y = 1.01, a x = 0.2475, a y = 0.2525, o/ x>y = 0.125, 
c = 0.4 and d = 0.01. 
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FIG. 12. Snapshots as in Fig. 11 at time t = 3200. 
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(d) 



FIG. 13. Snapshots as in Fig. 11 at time t = 17600. 
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FIG. 14. A snapshot at time t = 50500 of the field A x spontaneously generated from random 
initial conditions close to the trivial steady-state in the asymmetric case, a), b), c) and d) show 
the intensity, far field, real part and imaginary part of A x respectively. Parameters are as in figure 
11 except F = 1.06. 
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